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Abstract. - We investigate the finite time thermodynamics of a single-level fermion system 
interacting with a thermal reservoir through a tunneling junction. The optimal protocol to extract 
the maximum work from the system when moving the single energy level between an initial higher 
value and a final lower value in a finite time is calculated from a quantum master equation. The 
calculation also yields the optimal protocol to raise the energy level with the expenditure of the 
least amount of work on the system. The optimal protocol displays discontinuous jumps at the 
initial and final times. 



Introduction. — The search for the least work- 
intensive protocols for the extraction or insertion of en- 
ergy into or out of a thermal system has been a major 
research topic since the inception of the laws of thermody- 
namics. While the regime of quasi-static transformations, 
which is described by close-to-equilibrium thermodynam- 
ics, is well understood, many questions remain unsolved 
when dealing with problems far from equilibrium. A first 
question of particular interest deals with thermodynamic 
processes taking place in a finite time. This issue has 
been the object of detailed investigations in the context 
of finite time thermodynamics (FTT) [1]. Other devel- 
opments are related to recent progress in nanotechnology 
and cellular biology, where small systems far from equi- 
librium are subject to large thermal fluctuations. Devia- 
tions from average behavior, even rare events, play a sig- 
nificant role in their behavior. During the past decade, 
major progress has been achieved toward understanding 
and describing the role of fluctuations in such small non- 
equilibrium systems. The fluctuation theorem [2,3], the 
Jarzynski equality [4], Crook's theorem [5], and the for- 
mulation of stochastic thermodynamics [6] provide a novel 
framework to tackle the role of fluctuations in entropy pro- 
duction and dissipative work far from equilibrium. In ad- 



dition, exact expressions for the irreversible entropy pro- 
duction have also been derived [7-10]. A third frontline 
of research deals with quantum mechanical behavior in 
FTT. As the size of a system is reduced to the nanometer 
scale, quantum mechanical properties such as discreteness, 
quantum coherence, quantum statistics, and quantum cor- 
relations (entanglement) must be taken into account. We 
cite in particular the thermodynamics of quantum infor- 
mation processing [11,12], the related quantum heat en- 
gines [11, 13, 14], and quantum entanglement as a source 
of canonical typicality [15]. 

One of central questions addressed in FTT is to identify 
the optimal procedure to extract the greatest amount of 
work from a device operating under given constraints, or 
in reverse, to cause a device to operate under such con- 
straints with the minimum injection of work. According 
to the convention in which W is the work done on the 
system, maximum work extracted or minimum work in- 
jected both correspond to the minimum W. The question 
is thus that of identifying the protocol that involves the 
minimum amount of work done on the system. For exam- 
ple, Schmiedl and Seifert [16] considered the optimal pro- 
tocol to relocate a Brownian particle using a laser tweezer. 
They found that the optimal variation of the laser inten- 
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Fig. 1: (Color online) The model: a single-level quantum sys- 
tem interacts with a metallic thermal reservoir through a tun- 
neling junction. The level is initially at so and in thermal 
equilibrium. Work is extracted when the level is lowered to z\. 



sity which minimizes the work done on the system exhibits 
sudden jumps. Such singularities in the optimal proto- 
col may seem surprising, but in fact they turn out to be 
generic [17-20]. 

In the present letter we address a similar question for 
a simple quantum process. We consider a single level 
quantum system interacting with a heat bath. By rais- 
ing or lowering the energy of this level, we can inject work 
(W > 0) or extract work (W < 0). The time dependence 
of the protocol e(t) can be controlled externally. Our aim 
is to find an optimal protocol, one that minimizes the work 
done on the system, under the constraints of given initial 
and final values eo and e\, and a fixed total operation 
time r. We will specifically consider a quantum dot and 
a tunneling junction to a metal lead, the latter playing 
the role of a fermionic thermal reservoir. The detailed 
analysis of time-dependent phenomena in open quantum 
system is extremely complicated. In order to obtain ex- 
act analytical and numerical results, we restrict ourselves 
to a simple model based on a quantum master equation. 
We thus neglect quantum coherency and entanglement be- 
tween the system and the reservoir, but take into account 
the discreteness of the level and the proper Fermi-Dirac 
statistics. 

The model. — We consider a small quantum system 
consisting of a single level interacting with a metallic ther- 
mal reservoir through a tunneling junction, as illustrated 
in Fig. 1. A quantum dot with a tunneling junction to a 
lead provides a realization of such a system, provided the 
dot has a single energy level near the Fermi level of the 
lead, and no direct transitions between levels in the dot 
take place upon perturbation. We further assume that 
the electrons thermalize instantaneously upon tunneling 
to the reservoir, and that the latter remains in thermal 



equilibrium at a constant temperature T at all times. 

Upon varying the energy level or the chemical potential, 
a certain amount of (positive or negative) energy flows 
into the system as work done by an external agent. When 
the level is lowered (downward process) , work is extracted 
(W < 0). On the other hand, when the level is raised 
(upward process), work is injected into the system (W > 
0). The amount of work depends on the way the energy 
level is varied. Our goal is to find an optimal way of 
varying the energy level, the so-called optimal protocol, 
such that the maximum amount of work — W is extracted 
from the system or the minimum amount of work W is 
injected. 

We describe the time evolution of the quantum state 
using a master equation for the occupation probability 

p(*)> 

p(t) = -ui(t)p(t) + wa(t)[l - pit)], (1) 

where the are transition rates [21]. In the wide-band 
approximation, these rates are given by: 
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where C is a constant. Noting that raising the energy level 
is equivalent to lowering the chemical potential, we intro- 
duce an effective energy e(t) = s(t) — fJ,(t). By measuring 
time in units of C _1 , the master equation (1) thus reduces 
to the simple form 



p{t) = -p{t) + 



1 



1 



(3) 



We assume that the system is initially in thermal equilib- 
rium, 

*(°) - • (4) 
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Thermodynamic quantities. We next turn to a 
thermodynamic analysis of the model. We use the conven- 
tion that heat entering the system is (like work) positive. 
The internal energy of the system at a time t is 



where 



E(t) = U(t) - nN(t) = e(t)p(t), 



U(t)=e(t)p(t), N(t)=p(t). 



(5) 



(6) 



The rate of change in the internal energy, E, is the sum 
of two parts, namely a work flux VV and a heat flux Q, 



W — ep = ip — (ip 
Q = ep = ep - lip. 



(7a) 
(7b) 



Note that the particle exchange contributes to the heat 
flux [last term in Eq. (7b)]. When the energy level is below 
the Fermi level, the direction of heat flow is opposite to 
the direction of tunneling. 
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The net total work and net total heat during the process 
of duration r are obtained as functionals of the occupation 
probability, 



W[p(-)} = f e(t)p(t)dt 
Jo 

Q[p(-)} = f e(t)p(t)dt. 
Jo 



(8a) 
(8b) 



The resulting total net change in the internal energy is 
given by the First Law of thermodynamics, 

AE - p(r)ei - P (0)e = W\p(-)] + Q\p{-)] . (9) 

While work and heat depend on the path of p(t), AE 
depends only on the final probability p(r) and the given 
constraints e 0j e i andp(O). 

Minimizing work. 

General Approach. Our aim is to find an optimal pro- 
tocol e(t) which minimizes the work W. However, per- 
forming a variational analysis directly with respect to e(t) 
is complicated due to the expected discontinuities. In- 
stead, we optimize the work with respect to p(t), and 
identify the corresponding optimal e(t) from it. 

From the First Law of the thermodynamics, Eq. (9), we 
find work as a functional of p{t) , 



W[p(-)} = AE - Q[p(- 



(10) 



Since by definition p(t) is always differentiate, this ex- 
pression is well denned. 

In order to minimize work, we need to minimize AE 
and maximize Q simultaneously. However, from Eq. (9) 
we see that AE only depends on the final probability p(r) . 
Hence, we first identify the protocol leading to maximum 
heat Q for a given value of p(r). In a second step, we 
perform the optimization with respect to the final state 
p(r). To simplify notation, we assume in this section that 
energy is measured in units of kT. 

To find the protocol that maximizes the heat, we express 
e(t) in terms of p(t) and p(t) and rewrite Eq. (8b) as 

q\p(-)]= r c(p,p)dt, (a) 

Jo 



where 



£ = ln 



p(i). 



(12) 
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The extremum is found via the standard Euler-Lagrange 
method, leading, after integration, to 
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where K is the constant of integration. Before turning to 
the solution of this differential equation, we show that it 
implies a discontinuity in the protocol e{t). Eliminating p 



in Eq. (13) by using the master equation (3), the resulting 
quadratic equation for p{t) leads to the relation 
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1 ± VKe*W 



(14) 



If one determines the value of the integration constant K 
from the initial conditionp(O) assuming lim t _>o e(t) = e(0), 
this relation implies that K = 0, i.e., that p(t) is the 
equilibrium distribution associated with the instantaneous 
value of the energy. However, one expects that p{t) will 
deviate from thermal equilibrium except for an infinitely 
slow quasi-static process, so that in general K ^ 0. This 
apparent inconsistency indicates that lim t ^o e(0 7^ e (0). 
In other words, there must be a sudden jump from eo to 
e(0+). By comparing Eq. (4) to Eq. (14) at t = 0, we find 
the magnitude of the jump, 



s(0+) 



In 



1 ±2^ cosh 2 — 
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K COsh 2 j 
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Equation (14) also indicates that when K > there are 
two possibilities. The plus sign in ± leads to an occupa- 
tion probability p(t) that is larger than that of thermal 
equilibrium, and corresponds to the scenario of moving to 
a higher energy ei > eo- We refer to these as upward pro- 
cesses. For downward processes, the lower sign should be 
used. Henceforth it should thus be understood that the 
upper (lower) sign has to be considered when processes 
are upward (downward), respectively. 

Proceeding with the discussion of Eq. (13), we solve the 
quadratic equation for p, leading to 



K(l - 2p) T y/K 2 + AKp{\ - p) 
2(1 + if) • 



(16) 



This equation can be solved by separation of the variables 
t and p, leading to the following explicit result for the 
inverse function t as a function of p, 



t = F\p(t)]-F\p(0)}, 



(17) 



with 



F(p) 



= -ilnb(l- P )]T^arcsin(i=^ 

^arctanhPgr 1 ?! 4 ^"!^^!-^) 
2 \ 2(K- l)p(p- 1) + K J V ; 



While in general we will need to proceed with a numerical 
inversion for the resulting transcendental equation, an an- 
alytically tractable approximation will be discussed in the 
next section. Having thus obtained the optimal p(t) for a 
given K, we insert this expression in Eq. (8b) to obtain 
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the corresponding heat, 

rr rp(r) 

Q= / e(t)pdt = / e{p)dp 

JO Jp(0) 



cp(t) 



f 

= J dp In 

ip(0) 

= G\p(r)]-G\p(0)] 



K + 2p-2p 2 ± ^K 2 + AKp - AKp 2 



2p 2 



with rj(t) implicitly defined via Eq. (14). Next, we consider 
Eqs. (18) and (20) as functions of 77. Since r](t) -C 1 and 
e(t) <C 1 in the high temperature limit, we keep only lowest 
order terms. Noting from Eq. (13) that \[K is of the same 
order as e and 77, we find that Eqs. (25), (14), (18) and (20) 
simplify as follows: 



(19) 



where 



G{p) = pin 



K + 2p-2p 2 ± y/K 2 + AKp - AKp 2 



ln(l — p) ± v K arcsin 
2p - 2 



=F arctanh 



y/K 2 + AKp - AKp 




(20) 



Finally, we need to optimize the resulting work, given in 
Eq. (10), with respect to p(r), as explained earlier. Since 
p(r) is uniquely determined by K, it suffices to numerically 
optimize the expression with respect to K. 

We note some useful symmetries that arise from the 
fact that the optimal protocol to lower the energy level 
(with a resulting extraction of work from the system) is 
the mirror image in time of the optimal protocol to raise 
the energy level (associated with an injection of work into 
the system). Denoting the optimal protocol for upward 
processes from e to e\ > e by e-f (t) and the optimal 
protocol for downward processes from — e to — e\ by ei(t) 
and the associated occupation probabilities by p^ (t) and 
Pl(t) one easily finds that 



e T (i) + ei(t) = 0, 
p T (t) + Pi (*) = !■ 



(21) 
(22) 



These symmetries also imply symmetries for the minimum 
work and the associated heat, 



W^-W i 



Qi- 



(23) 
(24) 



These symmetries can be thought of as electron-hole sym- 
metries. They are particularly useful when the initial and 
final levels are symmetric with respect to the Fermi level, 
i.e., when e\ = — cq. 

The high temperature regime. The mathematical ex- 
pressions for the general case derived in the previous sub- 
section are rather complicated. However, the functions 
(18) and (20) simplify in the high-temperature limit, al- 
lowing us to find the optimal protocol and its properties 
in full analytical detail. 

First, we introduce an effective energy level r](t) defined 
by 

Pi*) = 377^ ( 25 ) 
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(26) 
(27) 
(28) 



T I^(t)-|„ 2 (t). 



Solving Eq. (17), we find rj(t) 
optimal protocol reads 



(29) 

e ± 2\[Kt and hence the 



e(t) = e a ±2VK(t + 1). 



The work is minimum for 



K= 1 -' 

A 



t + 2 

The optimal work and associated heat thus become 
(ei - e )[8 - 4e + r(4 - e 1 - e )] 



(30) 



(31) 



W 



Q = 



ir + 2) 



r{4 



8(r + 2) 



, (32a) 
(32b) 



From Eq. (30) we find that the initial and final energy 
jumps are given by ±f2\J~K. The size of the jumps in- 
creases as the deviation from the quasi-static limit (mea- 
sured by K) increases. In between the jumps, the op- 
timal protocol raises/lowers the level linearly with time 
(but we stress that this linear dependence only applies to 
the high temperature regime). Note that the results of 
the high temperature approximation satisfy the symme- 
try relations Eqs. (21)-(24). Here for the case with the 
symmetry eo = — ei there is no net heat flow. In this case 
all the work is converted into internal energy. 

Results. — In this section we present results for the 
optimal protocol obtained via numerical solution for a 
number of representative cases. First we consider the sit- 
uation where the energy level is raised from e = — lOfcT 
to ei = +10fcT during a total available time t = 10. 
This is the situation in which work is done on the quan- 
tum dot. The chosen parameter values guarantee that the 
initial energy level is well below, and the final level well 
above, the Fermi level. The top panel in Fig. 2 shows 
the optimal protocol. We also include two other protocols 
which do not have the optimal value of the initial jump, 
that is, ones corresponding to non-optimal values of K. 
The middle and bottom panels show the corresponding 
occupation probabilities and heat current. When the ini- 
tial jump is "too small" (dot-dashed lines), the level stays 
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Fig. 2: (Color online) Protocol (top), occupation probability 
(middle), and heat flux (bottom). The energy level is raised 
from e = -10 kT to ej = +10 kT during time r = 10. The 
red solid lines indicate optimal protocol with minimum work 
W = +10.87 kT. A protocol with a large initial jump (green 
dashed lines) and another protocol with a small initial jump 
(blue dot-dahed lines) result in higher work W = +16.79 kT 
and W = +11.16 kT, respectively. 



mostly below the Fermi level and thus the system receives 
heat from the reservoir. While this favors the reduction of 
work, there is not enough time for electrons, being below 
the Fermi most of the time, to tunnel, thus costing a large 
amount of work when raising the electrons during the fi- 
nal jump. On the other hand, when the initial jump is too 
large (dashed line), the electrons quickly escape and thus 
almost no work is required at the final jump. However, the 
rapid tunneling induces a large heat current to the reser- 
voir, increasing the work during the process. The optimal 
protocol (solid line) guarantees that the level is most likely 
empty before the final jump, with only a small outbound 
heat current. Note that the heat initially flows into the 
system, compensating for part of the outgoing heat after 
e(t) crosses the Fermi level. 

Figure 3 shows downward processes where the energy 
level is lowered from eo = +10 kT to ei = — 10 kT over 
the same period of time r = 10. In this scenario, work 
is extracted from the quantum dot. Three different cases 
parallel to the upward cases of Fig. 2 are illustrated. All 



15 

10 

5 



-5 

-10 

-15 
1 



a 0.5 




0.2 

3 

o 

^ -0.2 
-0.4 



W=-9.13 (optimal) 

W=-5.92 

W=-8.46 




2 4 6 8 10 
t 

Fig. 3: (Color online) Protocol (top), occupation probability 
(middle), and heat flux (bottom). The energy level is lowered 
from /3e = +10 kT to ei = -10 AT during time r = 10. The 
red solid lines indicate optimal protocol with minimum work 
W = — 9.13 AT. A protocol with a large initial jump (green 
dashed lines) and another protocol with a small initial jump 
(blue dot-dahed lines) result in higher work W = —8.58 kT 
and W — —5.92 kT, respectively. All three values of the work 
are related to the work in Fig. 2 through the electron-hole 
symmetry (23). 



the data shown in Fig. 3 confirm the electron-hole sym- 
metries surrounding Eqs. (21)-(24). 

Next, we reduce the time of operation. Figure 4 com- 
pares optimal protocols for slow (r = 10), intermediate 
(r = 1) and fast (r = 0.1) processes of work extraction 
from the quantum dot. As r decreases, the initial and final 
jumps become larger. When the electrons have almost no 
time to tunnel into the system, the optimal protocol be- 
comes nearly a step function, i.e., a jump-stay-jump pro- 
cess. 

Finally, Fig. 5 shows the optimal protocols for work ex- 
traction at various temperatures. As temperature is in- 
creased, the optimal protocol becomes symmetric with re- 
spect to the Fermi level, in good agreement with the high 
temperature approximation, Eq. (30). 

Discussion. — Discontinuities in the protocol that 
minimizes the work on a device operating under given 



p-5 



M. Esposito et al. 



1 
















T=0.1 






- > 


< 


T=l 
1=10 


































f 


1,1,1,1,1, 



0.2 0.4 0.6 0.8 1 



t 

Fig. 4: (Color online) The optimal protocols for three different 
processing times, t — 0.1 (solid line), r = 1 (dashed line), and 
t = 10 (dash-dotted line). The energy level is lowered from 
/3e = +10 kT to ei = -10 kT during time r. 

constraints may seem to be surprising. However, a simple 
phcnomcnological argument explains the initial and final 
jumps in the optimal protocols, starting from an analysis 
of the low temperature case. The key points to keep in 
mind arc: (i) moving an empty level requires no work, (ii) 
tunneling at the Fermi level carries no heat, and (iii) the 
direction of heat flow changes at the Fermi level. Con- 
sider an upward process at T = 0. Since tunneling is not 
possible below the Fermi level, the same amount of work 
is required to raise the level to the Fermi level regardless 
of the protocol. Hence, the instantaneous jump to the 
Fermi level is preferred since it leaves maximum time for 
electrons to subsequently tunnel out. After the jump, it 
is clear that the optimal protocol must keep the level in- 
hnitesimally above the Fermi level until the final time. In 
this way, heat transfer to the reservoir is avoided. Dur- 
ing this period no work is done, and the population in the 
energy level is reduced to close to zero without heat trans- 
fer. At the end, the level jumps up to its final value with 
almost no work. 

When the temperature is finite, tunneling is possible 
even below the Fermi level. Note that heat flux is inward 
into the system when the electrons tunnel out below the 
Fermi level, which helps reduce the work. On the other 
hand, the tunneling rate is small, which increases the work 
at a later time. The optimal protocol now includes an ini- 
tial jump to a level slightly below the Fermi level. At this 
stage, heat flows into the system. Next, the energy level 
moves slowly above the Fermi level. As electrons tunnel 
out, heat flows to the reservoir, which compensates the ini- 
tial heat gain. At the high temperature limit, the loss and 
gain of the heat are exactly balanced and no net heat flows 
to the reservoir during the process. All results shown in 
the previous section are consistent with this phenomeno- 
logical argument. 

We close with a critical discussion concerning the dis- 
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Fig. 5: (Color online) The optimal protocols for four different 
temperatures, T = 10 (green solid line), T = 1 (blue dashed 
line), T = 0.5 (black dash-double-dotted line), and T = 0.1 
(red dash-dotted line) . The energy level is lowered from eo = 1 
to £l = — 1, 

continuities in the optimal protocol. The jumps have been 
derived in the context of the master equation, which is 
valid only on a coarse grained time scale. Furthermore, 
we know from the time-energy uncertainty principle that 
an instantaneous jump in energy level would redistribute 
electrons over all energy levels including continuous states, 
i.e., the single level model cannot hold in this strict jump 
limit. Therefore, the discontinuities in the optimal pro- 
tocols identified here should be interpreted as rapid but 
continuous changes of the energy level. The correspond- 
ing typical time St must be much shorter than the tun- 
neling time (St <C 1). On the other hand, the single level 
model will remain valid only when St 3> 1/Aw, where Au 
is the gap between the energy levels. Even if the chemical 
potential is changed, we have assumed that the reservoir 
remains in thermal equilibrium. Hence, St must be longer 
than the relaxation time of the reservoir. These conditions 
can be satisfied if the tunneling rate is small, which is also 
a requirement for the validity of the master equation itself. 
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